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Single crystals of the organic charge-transfer salts a-(BEDT-TTF)2MHg(SCN)4 have been stud- 
ied using Hall-potential measurements (M =K) and magnetization experiments {M = K, Rb). The 
data show that two types of screening currents occur within the high-field, low-temperature CDWa, 
phases of these salts in response to time-dependent magnetic fields. The first, which gives rise 
to the induced Ifall potential, is a free current (jf^co)' present at the surface of the sample. The 
time constant for the decay of these currents is much longer than that expected from the sample 
resistivity. The second component of the current appears to be magnetic (jmag), in that it is a 
microscopic, quasi-orbital effect; it is evenly distributed within the bulk of the sample upon satu- 
ration. To explain these data, we propose a simple model invoking a new type of quantum fluid 
comprising a CDW coexisting with a two-dimensional Fermi-surface pocket which describes the two 
types of current. The model and data are able to account for the body of previous experimental 
data which had generated apparently contradictory interpretations in terms of the quantum Hall 
effect or superconductivity. 

PACS numbers: 71.45.Lr, 71.20.Ps, 71.18.+y 
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I. INTRODUCTION 

In their rudimentary form, charge-density waves 
(CDWs) consitute a simple one-dimensional spin singlet 
groundstate in which the collective mode consists of a 
charge density modulated at a characteristic wavevec- 
tor Q — 2fcF. Here, fcp is the Fermi wavevector of 
the Fermi-surface section responsible for CDW forma- 
tion Compared to superconductors, their behav- 
ior in a magnetic field is relatively easy to predict, pro- 
viding perhaps the simplest example of a singlet ground 
state reaching the Pauli limit |^, i^, i^, iTj . This limit is de- 
fined as the magnetic field required to cause the energy of 
the partially spin-polarized normal state to become lower 
than that of the condensate. 

In spite of the simplicity of CDW systems, the large 
values of the CDW condensation temperatures Tp, typi- 
cally of order 100 K HQ, make the Pauli limit inacces- 
sible in standard laboratory magnetic fields for the vast 
majority of CDW materials Our reason for study- 
ing the title compounds, Q!-(BEDT-TTF)2A/Hg(SCN)4, 
is that they lie at the far lower end of the spectrum of 
transition temperatures, with Tp « 8 K |^ for M = K 
and Tp Ki 12 K for M = Rb, making them perhaps the 
most marginal of CDW systems. Consequently, they are 
the only CDW compounds in which the gap is sufficiently 
low for the primary CDW phase, CDWq, to have been 
shown to be Pauli limited by a relatively modest mag- 
netic field ofjjpH « 23 T for M = K and hqH w 32 T 
for M = Rb nil mill Hi- 
lt must be stated, however, that the marginal CDW 
properties of a-(BEDT-TTF)2A/Hg(SCN)4 render it 
somewhat beyond the predictive power of the standard 
theory For instance, the low transition temperature 



results in a CDW with a weak charge modulation that 
is vulnerable to fluctuations and sample-dependent ef- 
fects 0. Furthermore, the presence of large sections of 
the Fermi surface left ungapped by the CDW order 
(see inset to Fig. ^ causes the energy associated with 
the coupling of the magnetic field to orbital degrees of 
freedom of the itinerant electrons to rival the condensa- 
tion energy. This leads to a modification of the quantum 
osciUatory effects j^ lT3IT5| and to the possibility of field- 
induced CDW sub-phases H, |3 ■ 

At magnetic fields above the Pauli paramagnetic limit, 
a-(BEDT-TTF)2MHg(SCN)4 then crosses over into a 
new regime in which the transition temparature becomes 
even lower 0, 0, falling to between 2 and 4 K: a 
schematic of the phase diagram is shown in Fig.^ Since 
the conventional CDW phase, CDWq, is no longer sta- 
ble at high magnetic fields, the spin-up and spin-down 
electrons must instead form independent charge and spin 
modulations with momentum vectors Qi and Qi (^|0], 
leading to a different type of CDW phase. 

Experimental studies of this high-magnetic-field phase, 
denoted CDW^,, lead one to question whether it is a CDW 
at all. While conventional wisdom has it that CDWs 
constitute a class of narrow-gap insulators P, 01 , exper- 
imental studies find behavior that is reminiscent either 
of the quantum Hall effectlll E El H 113 or super- 
conductivity m m in l24j. These experimental find- 
ings include a sharp drop in the electrical resistivity at 
low temperatures, persistent currents and unusual Hall- 
voltage phenomena. Problems arise because the quan- 
tum Hall effect '25^ and superconductivity "26] have fun- 
damentally different origins. If one attempts to catego- 
rize interpretations of experimental data in terms of the 
quantum Hall effect or superconductivity (23|, one 



obtains two sets of conclusions that are in contradiction. 

The present paper seeks to reconcile the body of ex- 
perimental evidence ElEHIllIiiililllllililin 

exploring the predicted properties of a new type of quan- 
tum fluid, consisting of a CDW coexisting with a 2D 
Fermi surface. A simple model shows that the exchange 
of quasiparticles between these subsystems results in the 
sample partially screening time- and spatially-varying 
electromagnetic fields via two types of current flow; this 
result is supported by the experimental data in this work, 
which show that the screening currents observed in ear- 
lier magnetic measurements 0, 0, [23, 0, IsSi comprise 
two parts. The first is weighted mostly towards the sam- 
ple edges and can be considered a free current jfroe, fur- 
nishing a Hall-potential gradient. The second contribu- 
tion is analogous to the "magnetic currents" jmag used 
to represent the effects of localised orbital moments 
in producing magnetization. These currents, which are 
essentially of infinite duration p^ . represent a local, mi- 
croscopic current flow within the bulk of the sample; they 
occur if the surface scrrening mechanism fails, causing the 
CDW, instead, to be coerced by an external change in 
magnetic field into a non-equilibrium state. In this case, 
bulk currents are supported by CDW pinning forces and 
are not free. 

The remainder of this paper is organized as follows. 
Experimental details are given in Section m whilst Sec- 
tion lllll describes the Hall potential measurements which 
characterize the surface current density jfree- The "mag- 
netic" current density jmag is studied using torque mag- 
netometry in Section IIVI Conclusions are given in Sec- 
tion El for ease of reference, the model describing the 




microscopic mechanism that leads to jftc 
scribed in the Appendix. 
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FIG. 1: Charge-density Wave formation in q-(BEDT- 
TTF)2MHg(SCN)4 salts. A schematic of the phase diagram 
as a function of reduced temperature T/Tp and reduced mag- 
netic field H/H]i (based on Ref. j23['l. The second order tran- 
sition temperature Tp at _ff = is approximately 8 K in the 
M salt and 12 K in the M = Rb salt. Likewise, the first 
order "kink" transition field hqH a.t T — (roughly corre- 
sponding to the Pauli paramagnetic limit), is 23 T for M =K 
and 32 T for M =Rb. 'CDWo' refers to the proposed conven- 
tional CDW phase while 'CDW^' refers to an incommensu- 
rate phase that comes into existence at high magnetic fields. 
The inset shows a simplified schematic of the Fermi-surface 
in the repeated Brillouin zone representation (after Ref. 27]) 
which consists of open one-dimensional (ID) sections which 
run along ky and a closed two-dimensional (2D) section. (An 
accurate representation of the actual Fermi surface topolgy is 
given m Ref. p.) CDW formation is thought to occur for a 
characteristic "nesting vector" Q 13| (with a component Q 
along kx), causing the ID section to become gapped. 



II. EXPERIMENTAL 

The free electrical current distribution, jftco, was de- 
termined using the variant of the Corbino geometry de- 
scribed in Ref. |2ll | in which pairs of small graphite-paint 
contacts are placed on the outer edge and upper surface 
of aa-(BEDT-TTF)2KHg(SCN)4 sample. Currents were 
induced in the sample by applying a small sinusoidal os- 
cillatory field of amplitude fioH and angular frequency lu 
superimposed on the quasistatic field iiqH provided by a 
33 T Bitter coil at NHMFL, Tallahassee. Alternatively, 
currents were induced by sweeping the quasistatic field 
at a rate fio{dH/dt). Provided that the sample dimen- 
sions are much smaller than its skin depth (so that it is 
entirely penetrated by the changing field), the resulting 
Hall potential between the edge of the sample, of area A, 
and its geometrical centre is j21| 

% = ^(^)a^o^ (swept field); (1) 



and 

Fh = — 1 — I ujfioH (oscillatory field), (2) 
An \ p|| J 

where p\\ is an appropriate average of pxx and pyy and 
standard symbols for the components of the resistivity 
tensor are used. 

Under the conditions of almost complete penetration 
by the oscillatory field, the phase of Vn in Eqn. |21 with 
respect to H is the same as that of a voltage induced in 
an open loop (i.e. 7r/2). However, if significant screening 
occurs, Vn will also contain a component in quadrature 
to this. Using standard complex notation, we write 

Vii = V^- iV^, (3) 

where and are the dissipative and reactive 
(quadrature) components of Vu respectively. In our ex- 
periment, phase-sensitive detection techniques allowed 
and to be simultaneously recorded for analy- 
sis purposes; the signal from a multiturn pick-up coil 
mounted close to the sample was used to define the real 
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phase. The approximate distribution of the current was 
further determined by placing pairs of contacts arranged 
between the outer and upper surfaces of the sample of 
a-(BEDT-TTF)2KHg(SCN)4 of volume « 1 mm^ shown 
in Fig. m While the finite size of 100 ± 50 fim of the 
electrical contacts applied using graphite paint limited 
the spatial resolution of the experiment, the contacts do 
enable one to distinguish between different models. Dur- 
ing these experiments, temperatures down to 0.5 K were 
achieved using a plastic '^He refrigerator. 

The orbital magnetic currents, jmag: were determined 
by measuring the magnetic torque of a different M = K 
sample to that shown in Fig. [5] as well as a M = Rb 
sample. In each case, 9, the angle between the magnetic 
field H and the normal to the sample's conducting planes 
was restricted to small angles {9 < 20°); this avoided 
complications due to magnetic torque interaction |27j| . 
The torque was measured capacitively by attaching the 
sample to a 5 /xm-thick phosphor-bronze cantilever that 
forms one of the plates of a capacitor. Electrical contacts 
applied to the sample confirmed that jmag continued to 
persist after the electrical currents dissipate. In these 
experiments, temperatures down to 50 mK were provided 
by a dilution refrigerator. 



III. HALL POTENTIAL MEASUREMENTS 

Fig. O shows typical measured Hall potentials and 
due to an oscillatory field; in this case fioH = 2.6 mT, 
iJ/2n = 409.6 Hz and the temperature of the q;-(BEDT- 
TTF)2KHg(SCN)4 sample {A « 1.1 mm^) was 0.5 K. 
The Vh data in Fig. |2| correspond to the contact labelled 
3 in Fig. 12 as the contact on the upper face of the sample 
was very near to its geometrical center, we expect that 
the Hall potential will be close to the maximum value 
possible. Note that there was a small field-independent 
inductive pick-up contribution Vp^ to due to the open- 
loop area (<1 mm^) between the contacts. By comparing 
measured values at several different fields, it is possi- 
ble to infer that Vpn ~ -4.7 fiY. 

In the discussion that follows, we first attempt 
to analyse the data of Fig |3| in terms of conven- 
tional bulk currents within the whole of the q;-(BEDT- 
TTF)2KHg(SCN)4 sample; the analysis uses conductiv- 
ity and resistivity tensor elements which do not vary spa- 
tially. We shall find that whilst this conventional analysis 
describes the behavior of V^, albeit with some remark- 
ably low values of it is unable to account for the 



variation of with H. 



A. Analysis in terms of bulk currents 

In a reasonably isotropic quasi-two-dimensional metal, 
the components of the resistivity (p) and conductivity 





< 1 mm^ 



FIG. 2: A photograph of the sample of a-(BEDT- 
TTF)2KHg(SCN)4Used in the measurements of the Hall po- 
tential, showing three contacts on the upper surface of the 
sample. Contacts on the side are not in focus. 
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FIG. 3: Examples of the real and imaginary Hall potentials 
and measured with the third pair of contacts on the sample 
shown in Fig. |5| An oscillatory applied field with ^oH « 
2.6 mT rms and u}/2n = 409.6 Hz was used; T ~ 0.5 K. 
Arrows indicate the axis corresponding to each data set. 



(cr) tensors are related by the expressions 

(Jxy _ CT|| 



p\\ 



Cm 



(4) 



Thus, the ratio Pxy/ P\\ in Eqn.|21is equivalent to axy/(^\\] 
it is therefore likely that the asymptotic variation of 
at low fields is due to the divergence of 

(y-xy = g2D/PQH, (5) 

where g2B is the charge density of the 2D holes, as H ^ 
0. 

As H increases, Pxyl P\\ becomes smaller, leaving Vpu as 
the dominant contribution to the measured Y^. However, 
after the "kink" transition at pi^H « 23 T, I^ undergoes 
a resurgence and the quadrature component becomes 
significant for the first time. The maxima in occur at 
integral Landau level filling factors v — F/p^H, where 
F is the magnetic quantum oscillation frequency; at such 



4 



fields, the chemical potential is situated in a Landau 
gap. The peak values V^^ « 25 ^J^^ (at ixqH « 25 T) 



correspond to Pxy/ P 
lows us to use Pxy ' 
Q2T)/e = 1.6 X 10^6 
pie resistivity of pu 



42: 



the fact that p^y S> p\\ al- 
i l/oxv = PoHj Q2T> (see Eqn.lSJ with 
m"^ |27| to extract a minimum sam- 
K, 2 Ax 10^^ ilm at the maxima close 



to /^o-ff = 30 T. Such resistivity values are characteristic 
of a good metal at room temperature and are not too dis- 
similar from the results of AC susceptibility experiments 
on a-(BEDT-TTF)2KHg(SCN)4 Hfl. 

To examine the apparent bulk resistivity components 
in more detail, we rearrange Eqs. 0]in terms of p\\/pxy 
and Uxy to yield 



Pxy 



1 



'xy 



and 



"! + ( — ) 



(6) 



(7) 



These equations can then be used along with Eqns. [3 
and|Slto convert the experimental values of to appar- 
ent bulk in-plane resistivity components. 

The results of this procedure are shown in Fig. ^ the 
linear increase of pxy with H and quadratic increase of p|| 
with H in at low magnetic fields show that the method 
gives a behavior in accord with expectations. At higher 
fields, the deduced values of pxy behave in a similar man- 
ner to earlier direct measurements of oscillations in the 
Hall resistivity However, it is interesting to note 

that the in-plane resistivity is lower in Fig. ^ than that 
obtained from four-terminal methods ,19J . Owing to non- 
uniform current flow associated with sample imperfec- 
tions, direct measurements of the in-plane resistivity are 
often contaminated by the inter-plane component p^z-, 
which is several orders larger in magnitude 
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FIG. 4: The bulk resistivity tensor elements p^y and py es- 
timated from Vh using Eqs. El and (T = 0.50 K). Arrows 
indicate the axis corresponding to each data set. 
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FIG. 5: An example of the Hall potential Vh measured for 
the third pair of contacts in a slowly varying magnetic field 
(sweep rate dpoH/dt = 0.5 T'^ T = 0.50 K). (a) shows 
raw data for a single sweep, with arrows indicating the sweep 
direction. No filtering has been applied. The background 
non oscillatory component could be caused by variations in 
the contact potentials with magnetic field, (b) The difference 
(Vh.up — VH,down)/2 betwceu rising- and falling-magnetic-field 
data averaged over 5 sweeps. 



The unusually large value of Pxy/ P\\ — 42 and con- 
sequent low value of py estimated from the oscillatory 

field measurements are confirmed by setting H — 
and by sweeping the quasistatic field. Data obtained 
using podH/dt = 0.5 Ts^^ are shown in Figure [S] for 
both rising and falling magnetic fields. Insertion of this 
rate of change of field into Eqn. ^ along with the value 
Pxy/p\\ ~ 42 obtained in the oscillatory-field experiments 
(see above), yields V-^ ~ 1.8 /iV at p^jH w 25 T, in 
good agreement with the data shown in Fig.|Sl However, 
the magnitude of the Hall potential continues to increase 
with increasing field, rising to a value of Vh ~ 8 /iV at 
around 32 T (Fig. [SJ. This yields Pxy/p\\ = 180, corre- 
sponding to p|| « 6 X 10~® rim, which is several times 
lower than the resistivity of room-temperature copper. 
Such values for these parameters would not be unusual 
for a two-dimensional electron gas exhibiting the quan- 
tum Hall effect 25]. However, they are unusual for a bulk 
organic metal far from the quantum limit; for example 
1/ w 21 at 32 T. 
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B. Analysis in terms of edge currents 

Thus far, a conventional treatment employing uniform 
conductivity tensor elements a\\ and a^y has given an 
apparently satisfactory explanation of the Hall poten- 
tial measured in swept fields and the real component 
of the Hall potential observed in oscillatory fields, 
albeit with extremely low values of the bulk resistiv- 
ity. However, this simple model fails on considering 
the imaginary component V{1 which appears within the 
high magnetic field regime in Fig. |3| Were q;-(BEDT- 
TTF)2KHg(SCN)4 a conventional metal, then in order 
to have a quadrature component of order 30 % of the 
in-phase component at a;/27r = 409.6 Hz, the skin-depth 
must be com para ble to or shorter than the effective ra- 
dius To = \/ A/tt ^ 590 /im of the sample Isij. If we 
estimate the skin-depth 5 — ^2p\^l [Iqlo from the above 
resistivity value (p|| ^ 2.4 x 10~^ fim), however, we ob- 
tain a value & = 3.9 mm that is significantly greater than 
the sample size. The presence of a large finite term, 
therefore, cannot be easily explained by a simple conduc- 
tivity tensor model. 

A further indication of the fact that the induced cur- 
rents do not obey a conventional bulk mechanism is given 
by using the contact arrangments shown in Fig. |21to ex- 
plore the dependence of and V{1 on the distance d of 
the contact on the upper surface of the sample from the 
edge. For a spatially homogeneous metallic system in a 
slowly time- varying magnetic field, it is trivial to show 
that the current density and concomitant Hall potential 
should vary quadratically; 



Vn{d) 



d{2ro - d) 



(8) 



where Vn.max is the total potential difference between the 
outer surface and geometric centre of the sample. This 
form occurs because the size of the area loop element sus- 
ceptible to induced voltages increases parabolically with 
the distance from the geometric centre of the sample. 
On comparing the detected Hall potentials for the three 
different contact arrangements shown in Fig. 12 we find 
that at ^qH « 25 T in Fig. varies more slowly 
with d than predicted by Eq. |S1 In the case of the reac- 
tive component, no discernable dependence of on d is 
seen. This type of behavior shows that the Hall poten- 
tial difference (and therefore also the current) is heavily 
weighted towards the sample edge, with the effect being 
particularly pronounced for the imaginary component. 

Currents carried at or close to the edge are expected to 
occur in quantum-Hall-effect systems for which cr|| — s- at 
integral Landau-level filling factors. These currents are 
either carried within edge channels ,25.] ■ giving rise to a 
kind of chiral Fermi liquid [33|, or are bulk states near 
the edge in which electrostatic forces locally deplete (or 
enhance) the charge carrier density |33| . Were a situation 
where a\\ reahsed in a-(BEDT-TTF)2KHg(SCN)4, 
perfect screening would give rise to a Hall potential that 
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FIG. 6: The Hall potential difference between the edge and a 
voltage probe situated a distance d inside the upper surface 
where d = 100, 310 and 400 /^m for contacts 1, 2 and 3 respec- 
tively, with T = 0.5 K, / = 409.6 Hz and fioH = 2.6 mT. (a) 
shows the normal Hall voltage V^h. The dotted line shows the 
voltage distribution expected according to Eq. |H| while the 
dashed line shows the voltage distribution expected for an 
exponential variation of the Hall potential, Vh oc exp(— d/A). 
(b) shows the reactive Hall voltage V^, with the solid line 
depicting Vh oc exp(— d/A) with A <C 100 /im. 
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FIG. 7: A comparison of the decay times tg^'^'^^ and to'^°'°*. 
All estimates were made for ^o-ff ~ 25 T and T = 0.5 K, 
with the data points corresponding to excitation frequencies 
lu/2tt = 51.25, 208.6, 409.6 and 8809 Hz and oscillating field 
amplitudes fioH = 6.51, 3.90, 2.62 and 3.75 mT. 



is entirely reactive, even in the limit as w — > 0. Ex- 
perimentally, |V[^'/^hI~ implying two things: first, 
CT|| ^ within the bulk, and second, the partial screen- 
ing that takes place has to occur within a region around 
the perimeter of the sample that is thinner than the bulk 
skin-depth. 

In the Appendix, we show that a system comprising a 
CDW coexisting with quasi-two-dimensional Landau lev- 
els can lead to very effective screening currents within a 
distance A — \/ m/2/ioeg2D of the sample surface. Al- 
though the finite size of the electrical contacts in Fig. [3 
(100 ± 20 /im in diameter), limits the resolution to which 
we can analyse the Hall voltage distribution in Fig. El 
we can, nevertheless, perform some simple consistency 
checks on this model. We propose that all of the cur- 
rent flows within a distance A of the surface of the sam- 
ple, and that this current flow encounters a characteris- 
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tic (very low) resistivity p\. The system is like an LR 
(inductance-resistance) circuit, with 



and 



L 



R 



I 



27rropA 



(9) 



(10) 



where I is the sample height, tq — yApK is its effective 
radius and a is a so-called Nagaoka par ameter, a numeri- 
cal factor '--^ 1 that depends on I/tq [SJI- For our sample, 
a ~ 0.5. In a magnetic field swept at a slow uniform 
rate, this treatment yields 

(c.f. Eqn. nj whereas for a sinusoidally-oscillating field 
we obtain 



roX f Pxy 



^(l-i-to). (12) 



PA / 1 + ^^t'o 



Here 



to = L/R = apQroX/2px 



(13) 



is a characteristic inductive decay time and Vnoe 
represents transient solutions due to rapid changes in 
dH/dt and/or R. 

We first compare Eqn. with the Hall potential 
of 8 /xV observed in Fig. [S] at pqH ~ 32 T, obtain- 
ing Pa/A w 23 ^ri, or Pa ~ 9 x 10^^^ fim, using 
A « 400 nm (see Appendix). Inserting this into the 
semiclassical Drude expression pA = 'm/eg^DT, we ob- 
tain a scattering time of r sa 50 ns, or, a mean free path 
A = upT ~ 5 mm (where we have obtained the Fermi 
velocity from vy = V2heF /m). The fact that A is com- 
parable to the sample perimeter of 4 mm could be con- 
sistent with ballistic transport. The ballistic transport 
regime, for which finite values of CT|| can no longer be lo- 
cally defined, is one of the essential preconditions of the 
edge- weighted current model (see Appendix). 

Two separate estimates of to can be made. First, we 
can rearrange Eqns. 1121 and US to yield to = aV^{l + 
LiPtf^) / PxyUjH . In the limit of low frequency, oj^tf^ ^ 1, 
so that this can be written 



PxyOjH 



(14) 



where the superscript denotes that this estimate of 
the characteristic time is derived from the dissipative (i.e. 
resistive) component of Vh- An alternative estimate can 
be made by comparing the resistive and reactive parts of 



decay 




(15) 



where the superscript '^'^'^'^^ denotes that this is analogous 
to the technique used to find the decay time of conven- 
tional LR circuits. 

Fig. H shows ^q"*'"'* and ig"'^^'*' deduced from experimen- 
tal data (po-ff ~ 25 T, temperature = 0.5 K). If the sam- 
ple surface layer functioned as a conventional (but high- 
conductivity) metal, we should expect the two times to 
be the same. However, the divergence between t'^'^'^^^ and 
^resist T^^-^ frequencies, or low values of Vjj, is unconven- 
tional. This suggests that the sample behaves resistively 
when attempting to drive a current through it, but once 
established, the currents have a tendency to last longer 
than expected for a conventional metal. Their lifetime 
^docay |-,ecomes particularly long at low values of the Hall 
potential Vyv observed in slowly- varying magnetic fields. 

In rapidly-varying magnetic fields, t'^^'^^^ converges 
with t]^^^^*- at T4j ~ 1 mV. This is similar to the satu- 
ration value of Vh observed in ms-duration pulsed mag- 
netic fields [21I I35I I . Attempts to explain this saturation 
in terms of the breakdown of the quantum Hall effect have 
been largely unsuccessful owing to the enormous magni- 
tude of the electric field that is required. An electric field 
that is concentrated towards the sample edge within a 
distance A (see Appendix), however, increases the proba- 
bility of Zener tunneling. Moreover, the increased size 
of the unit cell in the CDW phase compared to the 
normal metallic phase makes Zener tunneling across the 
CDW gap more likely than Landau level tunneling. In 
such a model for the magnitude of the threshold elec- 
tric field El = el/eaSF (where Eg = 2* is the CDW 
gap and sp is the Fermi energy), the lattice periodicity 
a = 2tt/Q « Tr/fcp occurs in the denominator. Assum- 
ing an exponential variation of the electric field where 
Vh = XEz, we obtain XEz = i'i'^X/nefiVF ~ 3 meV 
for 2^? = 1 meV, which is comparable to the saturation 
observed in pulsed magnetic field studies [21I Issj l . 

Should normal scattering processes be inhibited ac- 
cording to the screening model described in the Ap- 
pendix, Zener tunneling would provide a natural means 
for the decay of the currents, given the potentially large 
values of the electric field concentrated close to the sam- 
ple edge. 



IV. MAGNETIC TORQUE MEASUREMENTS 

The Hall potential associated with free currents jfree 
becomes difficult to detect for magnetic field sweep rates 
much less than 0.5 Ts~^. In contrast, magnetic currents 
jmag continue to be present for an arbitrarily long period 
of time after the magnetic field sweep is stopped [l^ . 
Fig.|Slshows examples of the magnetization of a-(BEDT- 
TTF)2MHg(SCN)4 (for M = K and Rb) estimated from 
the magnetic torque. Magnetic fields of poH ~ 23 T and 
Poi? « 32 T are required to access the CDW^, phase for 
M = K and Rb respectively [23|; in this phase, hysteresis 
similar to that originally obtained by Christ et al. p^ls^ 
is observed. 
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FIG. 8: The magnetization M of samples of the M — K 
{6 ~ 19.5°) and Rb {9 ~ 5°) salts measured on both rising 
and falling fields at T = 50 mK for M = K and T = 0.45 K 
for M = Rb made using the magnetic torque method. Arrows 
indicate the direction of sweep. For the M = Rb salt, at fields 
below approximately 29 T, the CDWo phase is stable. The re- 
gion between 29 T and 33 T appears to be disturbed since 
is different between rising and falling fields. At fields above 
approximately 33 T, the transition into the CDWa; phase is 
complete. The de Haas-van Alphen oscillations are the same 
for rising and falling fields, apart from an relative offset 2AM. 
For the M = K salt, the magnetic torque shows similar fea- 
tures to that measured in Reference The inset shows 
the magnetization measured after a field cool (as indicated in 
Fig. 1). The magnetization starts out at the value indicated 
by the large center-dot circle, after which the field is swept 
down to /io-ff < 29 T followed by an upsweep to ^nH > 30 T. 



The magnetization M of quasi-two-dimensional met- 
als such as a-(BEDT-TTF)2KHg(SCN)4 is largely di- 
rected normal to the conducting layers. A torque r = 
^^qMzH s\Tv29 per unit volume therefore occurs when 
there is a finite angle between M and and the applied 
field H • Sawtooth-like de Haas- van Alphen oscilla- 
tions (originating from the closed section of Fermi surface 
that survives the CDW formation p7|) account for the 
larger, reversible contribution to M. These oscillations 
occur in all metals with closed sections of Fermi surface 
at low temperatures and result from changes in the equi- 
librium populations of orbitally quantized Landau levels 
as H is swept [s^, IH, . The waveform of the oscil- 
lations in the M = K salt at high magnetic fields has 
received much theoretical attention and is well under- 
stood Hii^. 

By contrast, the contribution which we attribute to 
jmag is irreversible, and depends on how the sample is 
driven to a particular field and temperature. An exam- 
ple of this is shown in the inset to Fig. |S| if the sample 
is cooled to 0.5 K at fixed field, the magnetization at- 
tains the value surrounded by the circle. On sweeping 
the field down for the first time, the magnetization rises 
to meet the path which will be followed by subsequent 
downsweeps. Moreover, it proves to be impossible to go 
back to the encircled value of M by merely changing the 
field; data recorded on subsequent upsweeps fall below it, 




FIG. 9: Evidence for persistent currents in a-(BEDT- 
TTF)2KHg(SCN)4. (a) An example of a loop in the non- 
equilibrium component of the magnetization AM measured 
as the field is swept up to 29.33 T then momentarily down 
to 29.24 T before being resumed. Arrows indicate the change 
in the locus of AM versus H. The solid lines show the re- 
sults of fits of AM to Eq. 1171 (b) A model hysteresis loop 
calculated according to Eq. 1171 (i.e. the Bean model for a 
long cylinder with its axis parallel to H |40ll4H'). showing the 
theoretical saturation value of AM^ and the reversal state, 
(c) A series of AM versus H loops like that in (a) measured 
over an extended interval of field (top), and (bottom), a sim- 
ilar measurement but where the field interval over which the 
sweep direction is reversed is increased. Arrows indicate the 
appropriate axes, (d) Estimates of A^ from fits of Eg. 1171 to 
the many loops shown in (c) (full circles) together with those 
(large full squares) calculated according to Eq. 1321 (see also 
Ref. |2^). Estimates of the coercion field H* (x-symbols) are 
also shown. Arrows indicate the appropriate axes. 



whilst downsweep data lie above it (Fig. |S1 inset) . One 
can only recapture the field-cooled value by warming the 
sample well above the CDWa; phase and cooling again in 
the presence of a steady field. 

Thus, sweeping the field produces an irreversible 
contribution to the magnetization that manifests it- 
self as an offset 2AM between rising and falling field 
data in Fig Fig. shows AM for a a-(BEDT- 
TTF)2KHg(SCN)4 sample; values have been extracted 
by subtracting the reversible de Haas-van Alphen effect 
contribution (obtained by averaging full up and down 
sweeps) from the raw magnetization data to leave just 
the irreversible component. The magnetic field in Fig.|^ 
is first swept slowly up to 29.34 T, then down to 29.24 T, 
after which the up sweep is resumed. Each time the mag- 
netic field sweep direction is reversed, a finite interval in 
field AH > 2H* (where H* is the coercion field) is re- 
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quired in order for AM to reach a saturation value that 
opposes the direction of sweep. 

The Appendix describes how the CDW system pro- 
duces a diamagnetic contribution to the magnetization 
that acts to partiaUy screen the apphed magnetic field 
(see Eg. 1321 and the discussion following it). However, 
this contribution saturates when the force on the CDW, 
F, resulting from a spatially-varying free energy, exceeds 
the pinning force, Fp. Elastic energy considerations (see 
Appendix) show that the region in which |F| > |Fp| prop- 
agates inwards from the sample surface as the magnetic 
field changes, so that the volume- fraction of the sample 
able to contribute to the diamagnetic screening decreases. 
This is analogous to the scenario in the Bean model of 
type-II superconductors [29*, 40,'4l'|, yielding, under the 
approximation of a cylindrical sample, an irreversible sus- 
ceptibility 

Ax^^ = Axo/(Ai/) = Axo[l-^p. (16) 
and an irreversible magnetization 

rr* A rr 

AM=(Axo — )(2[1-— ]3-l). (17) 

Fig. ^jp shows the predicted behavior of the irreversible 
component of M according to Eq. El and Fig. shows 
a fit of the model to the experimental data. The fact that 
Eg. II 71 reproduces the experimental result indicates that 
the volume fraction of the sample f{AH) that is able to 
respond with an irreversible susceptibility Axo declines 
in a simple quadratic manner with AH (see Eq. , be- 
haviour tyracal of critical-state models such as that of 
Bean pfll l4ri l4ll | : here, f{AH) is the cross-sectional area 
of the remaining non-critical portion of the sample. 

As the magnetization varies throughout the sample, 
Maxwell's relation jmag = V x M [23 implies that a 
magnetic-field-induced critical state always involves cur- 
rents. However, although a critical-state model provides 
an excellent fit to the experimental data in Fig. it 
is important to note that there are two essential prop- 
erties of the currents giving rise to the loops in Fig. 
that distinguish them from those in superconductors. 
First, the magnitude of the non-equilibrium susceptibil- 
ity Axo ~ 3 X 10~^ (obtained from the fits) departs sig- 
nificantly from the ideal diamagnetic value Axo ~ — 1 
observed in superconductors. This implies that the mag- 
netic currents in a-(BEDT-TTF)2A/Hg(SCN)4 cannot 
be considered as screening currents in the same sense 
as those in superconductors. Because the magnetic cur- 
rents screen only a part of the magnetic field, the effective 
coercion field H* sa jcVo/Axo is rather large. Second, 
the value of A^o varies strongly as a function of H in 
Figs.Et andlHtl, a property that would be very difficult 
to explain in terms of superconductivity, but which is an 
intrinsic feature of the model of a CDW coexisting with 
well-defined Landau levels discussed in the Appendix. 

Fig. 1^ shows a series of loops in AM^ versus H like 
that in Fig. from which A^o is extracted at various 



values of H in Fig. |^. The magnitude of Axo is low 
and of order SxlO"** when /i is situated in the middle 
of a Landau level (at half- integral filling factors), but 
increases almost ten fold when /i is between Landau lev- 
els (at integral filling factors). The theoretical limit for 
Axo becomes equivalent to that of a superconductor (i.e. 
Axo = ~1) only in that case of an ideal sample at integral 
filling factors in which the quantum lifetime is infinite 
(see Eq. 1321 in the Appendix, remembering that demag- 
netising factors 123] must be taken into account once the 
susceptibility becomes very large). 

When pinning of the CDW occurs, the non-equilibrium 
susceptibility Axo provides a measure of how quickly the 
CDW departs from equilibrium as the magnetic field is 
swept. This occurs more rapidly at integral Landau level 
filling factors where the density of states is lowest, en- 
abling fi to jump quickly between Landau levels [29l |. 
Pinning of the CDW prevents the equilibrium redistri- 
bution of carriers between the bands, causing this jump 
to become much more abrupt. This can account the ori- 
gin of the strong variations in Axo in Fig. as shown 
by the model calculations (see Appendix). 



V. DISCUSSION AND CONCLUSION 

Before concluding this paper, it is intructive to 
compare effective ac susceptibilities of a-(BEDT- 
TTF)2KHg(SCN)4 determined using different methods 
in Fig. El In Fig. llOb . x' = ^Xo is obtained by fitting 
the Bean model (as described above- see Eq. I17II to many 
hysteresis loops over an extended region of magnetic field. 
The imaginary component x", which accounts for losses, 
is obtained from the loop areas. In Fig. llOb . the suscep- 
tibility is estimated using 

x'-^x"-^iVA-^V^). (18) 

The most clear aspect of Fig. El is that the magni- 
tudes of the real and imaginary susceptibilitiies due to 
free currents jfrco, estimated from the ac Hall potential, 
agree more closely with the high frequency ac susceptibil- 
ity measurements made in Ref. [s^ (shown in Fig. MUb ) 
than with the magnetic torque measurements. This in- 
dicates that the ac susceptibility made at frequencies 
uj/2'K 50 Hz consists mostly of free currents j tree, which 
are shown (above, and in the Appendix) to be confined 
to the sample edges. The data in Fig. |7| shows that jfico 
can persist for times as long as ^q'^'^^'^ 0.5 ms for suf- 
ficiently small values of the Hall potential Vh ^ 7 /iV 
in a sample of ^ 1 mm^ cross-section at ^qH « 25 T 
and T « 0.5 K. In spite of the fact that the correspond- 
ing susceptibility in Fig. [TUb is much larger than that 
for the torque measurements in Fig. llOb . the magneti- 
zation Mfrco = Vn/pxy ~ 7 Am~^ due to free currents 
is actually smaller. Free currents therefore have insuffi- 
cient duration to contribute a significant amount to the 
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FIG. 10: A comparison of ac susceptibilities (or equivalent ac 
susceptibilities) made on Q:-(BEDT-TTF)2KHg(SCN)4 using 
different methods. In (a) this is determined by fitting Eg. 1171 
to many hysteresis loops over an extended interval in mag- 
netic field. The imaginery component is determined from the 
area of the loop divided by its width. In (b), the effective sus- 
ceptibility is obtained by equating x' = 2nVi{ / pxyiJ,oH , where 
is the reactive component of the Hall potential. The imag- 
inary susceptibility is given by x' = 2-kV{i/ pxy^i-oH . Here, the 
susceptibility is estimated from the data in Fig. 3. (c) is the 
published ac susceptibility measured on another sample |30|. 



irreversible steady magnetic torque in Fig. |3 By con- 
trast, magnetic currents observed in the t orq ue experi- 
ments last for an indefinite period of time |l2l | . 

The combined induction of free and magnetic currents 
jfree + jmag Can therefore be understood as follows: for 
small changes in magnetic field hqAH <C 10 mT, the 
greatest contribution to the susceptibility initially origi- 
nates from free currents at the surface that attempt to 
screen changes in magnetic flux density from within the 
bulk. These currents quickly saturate, however, as dissi- 
pation sets in, possibly assisted by Zener tunneling. At 
that point, changes in magnetic flux density enter the 
bulk which cause the Landau level structure and chemi- 



cal potential to change, causing the CDW to depart from 
equilibrium, as described in the Appendix. The stored 
energy increases quadratically with AH, eventually caus- 
ing the CDW to collapse at the edges, initiating the crit- 
ical state. 

In conclusion. Hall potential and magnetic torque mea- 
surements on a-(BEDT-TTF)2MHg(SCN)4 (M = K, 
Rb) show that two types of screening currents occur 
within the high-fleld, low-temperature CDWj; phase in 
response to changing magnetic fields. The first, which 
gives rise to the induced Hall potential, is a free current 
(jfroc), weighted mostly towards the edge of the sample. 
The time constant for the decay of these currents is longer 
than that expected from the sample resistivity. The sec- 
ond component of the current appears to be magnetic 
(jmag), in that it is a microscopic, quasi-orbital bulk ef- 
fect; it is evenly distributed within the sample upon sat- 
uration. A simple model (Appendix), describing a new 
type of quantum fluid comprising a CDW coexisting with 
a two-dimensional Fermi-surface pocket, is able to ac- 
count for the origins of the currents. Taken together, 
these findings are able to reconcile the body of experi- 
mental evidence [IlllllllliillSilillllEl which 
had previously been interpreted in terms of the quantum 
Hall effect [23 or superconductivity p^ . 
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Appendix: Screening mechanisms in 
Q-(BEDT-TTF)2KHg(SCN)4 

We treat a simplified version of the a-(BEDT- 
TTF)2KHg(SCN)4 bandstructure 113, comprising a ID 
electron band and a two-dimensional (2D) hole band with 
dispersions Eid = hvF\kx — kp] and £20 — £f — Jt-"^ {\kx — 
fcxP + \ky — kY\'^)/2m respectively; here m is the effec- 
tive mass, wp is the Fermi velocity of the ID band, and 
kx and ky define the centre of the 2D hole pocket. Each 
of these bands intersects the Fermi energy (ep), giving 
rise to the simplified Fermi surface shown in the inset to 
Fig. 1. 

The simplest scenario to consider is that where only 
the ID Fermi-surface section is subject to CDW forma- 
tion, giving rise to a gap 25* in its density of electronic 
states, while the 2D hole section remains ungapped. Un- 
der equilibrium conditions, the average volume charge 
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densities giu and g2D associated with the ID and 2D 
Fermi-surface sections respectively would be subject to 
the conservation equation qid + ^2D + Qhg — 0, where gbg 
is the density of charge due to the ionic cores. However, 
below we consider slight local deviations ApiD and Ap2D 
in the charge density from the equilibrium values (i.e. the 
total local charge densities associated with the ID and 2D 
Fermi-surface sections become giu = ApiD + giu and 
g2D = Ap2D + g2D respectively); because of the presence 
of the CDW, the overall conservation equation need no 
longer hold locally. However, the charge densities must 
obey Poisson's equation 

-eVV(r) = Af?2D + ApiD, (19) 

where V is the electrostatic potential and e is the per- 
mittivity. We are at liberty to set the origin of poten- 
tial; we choose V = in the absence of spatial charge 
variations. Equation 1191 implies that the presence of 
local charge-density variations will lead to an in-plane, 
spatially-varying electric field E — —W. 

We now introduce the magnetic flux density Bq ap- 
plied along z (i.e. B = (0,0, B) = Bz, _L to the 2D 
{x,y) planes). This has two effects; first, the crossed E 
and B fields force the 2D hole wavefunction centres to 
drift at a velocity E x Hq/Bq. This leads to a spatially- 
varying in-plane current density, which, provided that 
the scattering rate is rather small (e.g. under conditions 
of ballistic or quasi-ballistic transport) , can be written as 
j « £>2dE X 'Bq/Bq-, note that the relationship becomes 
exact in the total absence of scattering. Integration of 
Maxwell's fourth equation (V x H = j shows that 
this current produces an additional contribution 

AB = -^iog2DiV/Bo)z (20) 

to the magnetic fiux density parallel to z, which then be- 
comes B = Bo -I- AB. A second effect of B is to produce 
Landau quantisation of the 2D holes; this is usually dealt 
with using the Landau gauge A = {0,Bx,0) Af- 
ter some manipulation |42j , the Schrodinger equation for 
the in-plane wavefunctions ip and eigenenergies E may be 

written as (^—j^-^ + ^muj'^{x — Xq^^ tp = Eijj. Here, 
the effect of the crystalline potential has been taken into 
account by the effective mass m; the cyclotron frequency 
is LOc — qB/m, with q the charge; xq represents the wave- 
function guiding-centre coordinate. 

The in-plane electric field also modifies the Landau- 
level energies [sj. To quantify this, we consider the case 
V(r) = V{x), i.e. the potential varies only in the x 
direction (variation in an arbitrary direction in the x, y 
plane is reintroduced below). If V(x) varies slowly over 
lengthscales ~ xo, it can be expanded about xq in a 
Taylor's series V{x) ~ V{xo) + {x — xo)V'{x — xq) + 
^{x — xq)^V" {xq) + , where the primes indicate dif- 
ferentiation with respect to x. After a little algebra, the 
Schrodinger equation becomes 

-i^P^ + ilM^! + -V"{xo)Kx~x^r)^ = E^, (21) 
Zm ax^ 2 m 



where xi is a constant (absorbing terms in xq, ^(a^o) E^nd 
V'ixo)). 

The chief effect of V{x) comes from V"{xo); this gives 
a modified Landau- level spectrum, {n + ■^ jhui, with n an 
integer and 

(22) 

where Wco = qBo/m and only terms to leading order in 
AB and I^" are retained in the right-hand bracket. Rein- 
troducing a potential variation in an arbitrary intraplane 
direction merely changes V" in Eqn. |221to W'^V. 

In the absence of potential variations, the number of 
states per unit volume per Landau level is D = Dq = 
rnwco /7r?ic, where c is the layer separation in the z direc- 
tion; the introduction of varying V (and consequent AB) 
causes uJqO to change to uj (Eqn. I22() , modifying this to 
D ~ muj /-kTic = D{) + AZ?, where 

AD D (^^ + ^ (23) 

\ Bq 2uJcaBoJ 

The change in Landau-level degeneracy results in a 
change to the local 2D hole density, 

Aei2D = i^qAD - vPqAD - f3qq2BA^i. (24) 

It is worth taking time to understand this equation fully, 
since it is the key to understanding the screening effects 
that are the point of this paper. The first term of the 
right-hand side results from the change in degeneracy of 
all of the occupied Landau levels; positive V^F (r) or AB 
gives a greater Landau-level degeneracy fEan. pro- 
ducing an increase in the number of holes. The second 
term results from the energy shift of the I'th Landau level 
closest to the chemical potential /i; in the high-field limit 
of well-resolved levels, the contribution from the tails of 
the Landau levels remote from /i can be ignored. As 
an increase in V^y(r) causes uj to increase (Eqn. I22|l . 
this results in a slight depopulation of the Landau levels 
closest to ^, acting to reduce p2D- The dimensionless 
factor P therefore depends critically on the position of fi 
amongst the Landau levels; i.e. it oscillates as a function 
of B. If /X is in the middle of the highest occupied Landau 
level (at half- integral Landau- level filling, v = F/ B), the 
density of states is large, so that the shift in energy of 
the Landau level has a relatively large effect on p2d; in 
this situation, /3haif ~ 2ujct/'it ^ 1. On the other hand, 
if /I is directly between two Landau levels (at integral 
Landau-level filling v = F/B) then the density of states 
is small and /3int ~ A/ttlOcT <^ 1. This tends to zero 
in the case of a sample of pristine purity for which the 
quantum lifetime r — > oo |33|. The final term accounts 
for possible spatial variations A^ in the chemical poten- 
tial which are sustained by the CDW (see below); here, 
^2D = g^DTTi / TiF q^ is the mean value of the 2D density 
of states. 
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A shift in the chemical potential also affects the ID 
carrier density; 



AgiD = -eginAfi, 



(25) 



where ijid is the ID density of states. We can now re- 
formulate Eqn.Elby substituting from Eqns. [1111101 ESI 
and 1251 to yield 



2w, 



-V^V + AB 



cO 



f3 + r] 
hev \ l- (3 



A/i. 



(26) 



We have retained terms up to first order in small quan- 
tities, omitted eV^V from Eqn. ^| because it is ~ 10"* 
times smaller than the other term in V'^V and substi- 
tuted q — +e for the hole charge; here, -q — (7id/.92D- 



2D hole Landau levels 




1 D electron dispersion 




FIG. 11: A model depicting the affect of a change A/^ oc 
AA'^iD in the chemical potential /x on the 2D hole Landau 
levels and the ID electron dispersions (in this case at half- 
integer filling), giving rise to a direct relationship between 
AA''iD, A^ and AQ and consequently between AQ and H. 

As long as the ID band remains metallic, A/i is un- 
constrained. Scattering facilitates the transfer of current 
between positive and negative variations in the charge 
density, causing them to dissipate rapidly. The stable 
result is when V'^V — 0, causing Eqn. |2Slto become a 
simple proportionality between A/i and A_B; i.e. a varia- 
tion in the chemical potential due to a change in magnetic 
field 

By forming a CDW on the ID bands, the system can 
benefit in two ways. First, by the opening of a CDW 
gap 2^*, the system has the potential to become partly 
gapped at /i. Second, by forming a quantum-coherent 
state, the ID band resists the spontaneous exchange of 
charge between the two bands. As with a superconduc- 
tor, the quantum state of a CDW is locally defined by a 
phase (/), which defines the phase of the charge modula- 
tion that oscillates in one direction [x) on a lengthscale 
27r/|Q| much shorter than the cyclotron length. Any shift 
in the chemical potential in Eqn. will correspond to a 
gradient in phase (see Fig. ITTl) ; 



A/x = 



Tivy d4> 



(27) 



The proportionality between d4>/dx and A/i leads to two 
qualitatively different solutions to Eqn. 1261 



Within the bulk of the sample 

In the bulk of the sample, the transport is dissipative 
(V^F 0) and the CDW is able to adjust its phase 
within certain limits to accommodate currents caused by 
a non- uniform potential. Eq. 1261 therefore becomes 



A/i _ heiy / 1 - /3 
AB ~ m \P + 7-i 



(28) 



Thus, there is a direct proportionality between local vari- 
ations in the magnetic field and A/i. As noted above 
(Eq. 123), a change A/i results in a shift in CDW phase, 
and, via 1231 



dAQ 
djj, 



,ffiDe 
IpidI 



(29) 



to local adjustments of Q away from its equilibrium value, 
Q (see Fig. Ill|l . These changes result in a concomitant 
shift in the free energy $ of the system, which, following 
Ref. |29|, can be written as 



-5iD 



^2 



A/i^ + 5gmAti^ + $2d(£'id). (30) 



Here, $2D is the free energy of the Landau-level system 
(which depends at fixed field on the number of quasi- 
particles in the 2D hole band, and hence on gm) and 5 
is a dimensionless parameter depending on the position 
of the chemical potential amongst the Landau-levels; it 
takes the limiting values [2^ 



1 



1 



27J^U;,T 2(77fWeT) 



(31) 



for integer Landau-level filling factors and (5 = ^ for half- 
integer filling factors. 

The effects described by Eqs. I28I30I result in a rela- 
tionship between $ and AB, implying that there is a 
contribution Ax from the CDW to the overall suscep- 
tibility X, where Ax = — /io(i9^<&/i9(Ai3)^). Using the 
identity {d'i>/d{AB)) = (9$/a(A/i)) x (a(A/i)/a(AB)) 
in conjunction with Eas. and IHUI vields ^44] 



Ax = -2<?iD/io(l + '5) 



hev 




')! 


m 







(32) 



Thus, the non-equilibrium (i.e. spatially-varying) state 
of the CDW thus provides a diamagnetic contribution 
to the overall susceptibility, which resists changes in the 
magnetic field in the bulk of the sample; note that the 
presence of the factors (3 and 5 indicate that the effec- 
tiveness of this screening will oscillate as a function of 
Landau-level filling factor u. However, as we discuss in 
the next paragraph, Eq.|221represents the response of the 
sample only to infinitessimal changes in magnetic field; 
the build up (and release) of elastic energy as the field 
develops further leads to an irreversible change in the 
overall sample susceptibility which eventually saturates. 
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The spatially-varying free energy associated with the 
mechanism causing Eq. corresponds to a force F — 
— V$; the non-equilibrium excited state of the CDW will 
only persist if F remains less than the CDW's pinning 
force Fp This produces two restrictions; first, be- 

cause of the differential relationship between F and A/i, 
Afi must tend to zero at the sample surface, or else 
F would be singular (and hence exceed Fp). Second, 
once F develops with changing field to such an extent 
that the CDW depins, the non-equilibrium build up of 
Afi must stop, preventing any further contribution from 
Ax- These combine to ensure that the region in which 
|F| > |Fp| propagates inwards from the sample surface as 
the magnetic field changes, so that the volume-fraction 
of the sample able to contribute Ax (Eq. I32|l to the 
overall volume-averaged sample susceptibility x becomes 
progressively smaller. This is very similar to the Bean 
critical-state model [23^3,^3; indeed, a common fea- 
ture of all critical-state models is that motion initiates 
at the surface because it is there that the potential en- 
ergy can be dissipated with the least amount of action. 
The CDW begins to slide at the sample surface so that 
the restoring force does not exceed the CDW depinning 
force. Approximately linear gradients in AgiD and AQ 
occur because the work done by F is a linear function 
of the departure AQ of the CDW from equilibrium ^2^. 
Currents result from changes in the charge density g2B 
associated with the two-dimensional Fermi surface sec- 
tion as the quasiparticles redistribute themselves so as 
to screen the changes in the open Fermi surface section 
(see Eq. I19II . The variation in the orbital magnetization 
of the closed Fermi-surface pocket AM with Ag2ri then 
leads to linear gradients in AM. Maxwell's equations im- 
ply that a linear gradient in AM, sustained by pinning, 
corresponds to a critical current density jc = V x (AM), 
in the region where the CDW is able to shde. 

Therefore, although the mechanisms are very different, 
there are close analogies between the situation here and 
fiux-pinning in superconductors. It is therefore not sur- 
prising that in Section llVl we shall see that a simple Bean 
model is able to fit the variation of sample susceptibility 



with applied magnetic field. 



Close to the sample surface 

We have already seen that Afj. — > as one approaches 
the sample surface because the elastic energy cost of shift- 
ing the phase would otherwise become too great. Hence, 
Eq. assumes the form 

X^V^V-V^O with A = [m/(2^oe^2D)]^ (33) 

where we have substituted for AB using Eq. [201 Eq. ISSI 

defines a penetration depth A analogous to the Lon- 
don penetration depth in superconductivity; it implies 
that deviations V from the equilibrium value of the elec- 
trostatic potential will be screened from the bulk of 
the sample. Owing to the proportionality between V 
and AB in Eqn. 1201 spatial variations in the magnetic 
field will also be screened. Substituting values for a- 
(BEDT-TTF)2KHg(SCN)4 taken from Ref. ^ yields 
A w 400 nm. 

The existence of a finite V'^V close to the surface im- 
plies quasi-ballistic transport. By carrying the current 
close to the edges, the total system can save energy 
(in this case, elastic energy due to the non-equilibrium 
arrangement of the CDW within the bulk of the sam- 
ple), and any process that saves energy also protects 
the quasipartices from ordinary scattering events that 
would otherwise give rise to dissipative bulk currents. 
This type of explanation for the variation in the cur- 
rent density, electrostatic potential, electric field and per- 
turbed charge carrier density follows on from a model 
originally proposed by MacDonald et al. [s^ for two- 
dimensional electron gas systems. Whereas for a sin- 
gle two-dimensional layer an exact solution can only be 
obtained numerically, in bulk crystalline systems like a- 
(BEDT-TTF)2KHg(SCN)4, the above treatment shows 
that the variation in all of these quantities becomes a 
simple exponential function. 
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